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OPTIMAL ARBITRAGE UNDER MODEL UNCERTAINTY 

By Daniel Fernholz and Ioannis Karatzas^'^ 
Daniel Fernholz LLC and INTECH Investment Management 

Dedicated to Professor Mark H. A . Davis on the occasion of his 65th 

birthday 

In an equity market model with "Knightian" uncertainty regard- 
ing the relative risk and covariance structure of its assets, we charac- 
terize in several ways the highest return relative to the market that 
can be achieved using nonanticipative investment rules over a given 
time horizon, and under any admissible configuration of model pa- 
rameters that might materialize. One characterization is in terms of 
the smallest positive supersolution to a fully nonlinear parabolic par- 
tial differential equation of the Hamilton-Jacobi-Bellman type. Un- 
der appropriate conditions, this smallest supersolution is the value 
function of an associated stochastic control problem, namely, the 
maximal probability with which an auxiliary multidimensional dif- 
fusion process, controlled in a manner which affects both its drift 
and covariance structures, stays in the interior of the positive orthant 
through the end of the time-horizon. This value function is also char- 
acterized in terms of a stochastic game, and can be used to generate 
an investment rule that realizes such best possible outperformance of 
the market. 

1. Introduction. Consider an equity market with asset capitalizations 
X(t) = . . . ,Xn{t))' G (0,oo)" at time t G [0,oo), and with covariance 

and relative risk rates a(t, X) = {aij{t^ ^)}i<i,j<n and -dit, X) = {"diit^ X), . . . , 
'dn{t,X))\ respectively. At any given time t, these rates are nonanticipative 
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Junctionals of past-and-present capitalizations X{s), < s <t; they are not 
specified with precision but are, rather, subject to "Knightian uncertainty." 
To wit, for a given cohection 

(1.1) 1K = {/C(y)}yee„, 6„ := [0, oo)" \ {0} 

of nonempty compact and convex subsets on M" x where S" is the space 
of real, symmetric, positive definite (n x n) matrices, and is the origin 
in M^, they are subject to the constraint 

(1.2) (?9(t, X),a{t, X)) G JC{X{t)) for all t G [0, oo). 

In other words, the pair (??,a) must take values at time t inside the com- 
pact, convex set IC{X{t)) which is determined by the current location of 
the asset capitalization process; but within this range, the actual value 
{'&{t,X),a{t,X)) is allowed to depend on past capitalizations as well. [To put 
it a little differently: the constraint (1.2) is not necessarily "Markovian," as 
long as the sets in (1.1) are not singletons.] 

Under these circumstances, what is the highest return on investment rel- 
ative to the market that can be achieved using nonanticipative investment 
rules, and with probability one under all possible market model configu- 
rations that satisfy the constraints of (1.2)? What are the weights in the 
various assets of an investment rule that accomplishes this? 

Answers: Subject to appropriate conditions, 1/U(T,X{0)) and 

(1.3) 

i = l,...,n,0<t<r, 

respectively. Here the function U :[0,oo) x (0,oo)" — ?• (0,1] is the smallest 
nonnegative solution, in the class C^'^, of the fully nonlinear parabolic partial 
differential inequality 

dU 

(1.4) —iT,z)>CU{T,z), (t,z)g(0,oo) X (0,oo)" 
subject to the initial condition U{0, •) = 1, with 

(1.5) £/(z)= sup x:e^.^.«./^4/(z)+ , )• 

We use in (1.3), (1.5) and throughout this paper, the notation Dif = df/dxi, 
Dfjf = d^f/dxi dxj, Df = {Dif, . . . , DnfY, D^f = {Dy}i<i,,<n and define 

(1.6) ^(y) := {a G S" : (9, a) G /C(y), for some 9 G M"}, y G ©„. 



We call the function U{-, ■) the arbitrage function, as U (T, x)(a;i H \-Xn) 

gives the smallest initial capital starting with which an investor, who uses 
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nonanticipative investment rules, can match or outperform the market port- 
foho by time t = T, if the initial configuration of asset capitalizations is 
X{0) = X = (xi, . . . , Xn)' S (0, cxd)" at t = 0, and does so with probability one 
under any "admissible" market configuration that might materialize. It is 
perhaps worth noting that this function U{-,-) is characterized almost en- 
tirely in terms of the prevalent covariance structure a. The relative risk en- 
ters only indirectly, namely, in determining the family of sets (1.6) which are 
admissible for the covariance structure. Put a bit differently, the only role 
plays is to ensure the asset capitalization process X(-) = (Xi(-), . . . ,X„(-))' 
takes values in (0,oo)"'. 

Under additional regularity conditions, U (T, x) is the value of a stochastic 
control problem: the maximal probability that the diffusion process 2)(-) = 
(Yi(-), . . -M-))' with initial configuration 2)(0) = X(0) = x e (0,oo)", val- 
ues in the punctured nonnegative orthant S„ of (1.1), infinitesimal generator 



and controlled through the choice of covariance function a : [0, oo) x — )• 
which satisfies a(t,y) G ^(y) for all (i,y) G [0,oo) x S„, does not hit the 
boundary of the orthant [0, oo)" by time t = T. Under appropriate conditions 
the function U{-,-) satisfies then, in the notation of (1.5), the Hamilton- 
Jacohi-Bellman (HJB) equation 



Relation to extant work: Stochastic control problems of the "maximal prob- 
ability of containment" type were apparently pioneered by Van Mellaert 
and Dorato (1972) (see also Fleming and Rishel (1975), pages 157-158). 
The "Knightian uncertainty" constraint imposed in (1.2) is very similar to 
the formulation of stochastic control and stochastic game problems for one- 
dimensional diffusions pioneered by William Sudderth, that Sudderth and 
his collaborators developed in a series of articles that includes Pestien and 
Sudderth (1985), Heath et al. (1987), Orey, Pestien and Sudderth (1987), 
Sudderth and Weerasinghe (1989); indeed, the developments in Sections 6-8 
of our paper can be construed as a multidimensional extension of the Sud- 
derth approach. 

We rely strongly on Krylov's (1989, 2002) work, which studies solutions of 
stochastic differential equations with constraints on the drift and diffusion 
coefficients in terms of "supermartingale problems" and characterizes sets 
of stochastic integrals via appropriate supermartingales. 

The approach we adopt has a lot in common with the effort, started in 
the mid-1990s, to understand option pricing and hedging in the presence 
of uncertainty about the underlying volatility structure of assets. We have 




(1.7) 



{d/dT)U{T,z) = CU{T,z) 



on (0,oo) X (0,oo)". 
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been influenced by this strand of work, particularly by the papers of Lyons 
(1995), Romagnoli and Vargiolu (2000), Gozzi and Vargiolu (2002), Vargiolu 
(2001), Talay and Zheng (2002); other important papers include Avellaneda, 
Levy and Paras (1995), El Karoui, Jeanblanc-Picque and Shreve (1998), 
Cvitanic, Pham and Touzi (1999), Frey (2000), Ekstrom and Tysk (2004), 
Meyer (2006), Denis and Martini (2006), whereas the recent preprints by 
Soner, Touzi and Zhang (2010a, 2010b) contain very relevant results. Sim- 
ilar in this spirit is the strand of work by Shige Peng and his collabora- 
tors, surveyed in Peng (2010), regarding the so-called "G-Brownian mo- 
tion" which exhibits volatility uncertainty [see also Vorbrink (2010), as well 
as Nutz (2010) for extensions to settings where the range of uncertainty 
is stochastic]. Whereas in both these strands the relevant fully nonlinear 
parabolic- type partial differential (so-called "Black-Scholes-Barenblatt" ) 
equation has a typically unique solution, here the main interest arises from 
lack of uniqueness on the part of the rather similar, fully nonlinear equa- 
tion (1.7). 

Let us mention that optimization problems in stochastic control, mathe- 
matical economics and finance that involve model uncertainty have also been 
treated by other authors, among them Gilboa and Schmeidler (1989), Gun- 
del (2005), Shied and Wu (2005), Karatzas and Zamfirescu (2005), Follmer 
and Gundel (2006), Schied (2007), Riedel (2009), Bayraktar and Yao (2011), 
Bayraktar, Karatzas and Yao (2011) and Kardaras and Robertson (2011) 
[see also the survey by Follmer, Schied and Weber (2009)]. 

Preview: Sections 2 and 3 set up the model for an equity market with 
Knightian model uncertainty regarding its volatility and market-price-of- 
risk characteristics, and for investment rules in its context. Section 4 in- 
troduces the notion of optimal arbitrage in this context, whereas Section 5 
discusses the relevance of the fully nonlinear parabolic partial differential 
inequality of HJB type (1.4), (1.5) in characterizing the arbitrage function 
and in finding an investment rule that realizes the best outperformance of 
the market portfolio. Section 6 presents a verification-type result for this 
equation. Sections 7 and 8 make the connection with the stochastic control 
problem of maximizing the probability of containment for an auxiliary Ito 
process, controlled in a nonanticipative way and in a manner that affects 
both its drift and dispersion characteristics. Finally, Section 9 develops yet 
another characterization of the arbitrage function, this time as the min-max 
value of a zero-sum stochastic game; the investment rule that realizes the 
best outperformance of (optimal arbitrage with respect to) the market, is 
now seen as the investor's best response to a "least favorable" market model 
configuration. 

2. Equity market with Knightian model uncertainty. We shall fix through- 
out a canonical, filtered measurable space (17,7-"), F = {^(t)}o<t<oo and as- 
sume that Q contains the space 2IJ = C([0, oo); (0, oo)") of all continuous 
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functions tr : [0, oo) — t- (0, oo)". We shall specify this canonical space in more 
detail in Section 7 below, when such detail becomes necessary. 

On this space, we shall consider a vector of continuous, adapted processes 
X(-) = (Xi(-), . . . ,Xn{-)y with values in (0,oo)"'; its components will repre- 
sent stock capitalizations in an equity market with n assets, and thus the 
total market capitalization will be the sum 

(2.1) X{t):=Xi{t) + --- + Xnit), 0<t<oo. 

We shall also fix throughout a collection K = {/C(y)}yge„ of nonempty, 
compact and convex subsets of M" x S" as in (1.1). 

We shall consider ]R'"-valued functionals •) = (??i(-, ■),■■■ ,'dn{', ■))' 
S"-valued functional a(-,-) = (aij(-, •))i<ij<n; of them defined on [0,oo) x 
0, and progressively measurable [see Karatzas and Shreve (1991), Defini- 
tion 3.5.15]. We shall assume that, for every continuous function ro : [0, oo) — >• 
(0, oo)" and T £ (0, oo), these functionals satisfy the constraint and integra- 
bility conditions, respectively, 

(i9(r,tt)),a(r,tt))) G/c(tt)(r)), 

(2-2) 

/ (||i9(t,m)f + Tr(a(t,w)))dt<oo. 
Jo 

We shall also consider (n x n) -matrix- valued functionals 
where 

a{t, ro) = \/a(t, id) is a square root of a{t, It)) : a{t, Id) = a{t, to)a'{t, ro). 

(2.3) 

2.1. Admissible systems. For a given collection of sets IK as in (1.1) and 
a fixed initial configuration x = {xi, . . . , Xn)' G (0, oo)" of asset capitaliza- 
tions, we shall call admissible system a collection Ai consisting of the un- 
derlying filtered space F = {J^{t)}o<t<oo, of a probability measure P 
on it, and of a pair of processes (X(-), W{-)), with W{-) = {Wi{-), . . . , W„(-))' 
an n-dimensional F-Brownian motion and X(-) = (Xi(-), . . . , X„(-))' taking 
values in (0,oo)". These processes have the dynamics 

n 

(2.4) dXi{t)=X,{t)^aUt,X)[dW,{t) + Mt,^)dt], XiiO) = x^>0 

for some progressively measurable functionals '&{■,■) and which sat- 

isfy (2.2) and (2.3) above. 

We shall think of this admissible system as a model subject to "Knigh- 
tian" uncertainty; this is expressed by the requirement (i?(t, jC), a(t, X)) G 
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/C(X(i)) in (2.2), (1.1) about the market price of risk and the covariance 
structure of the asset capitahzation vector process X(-). In order not to lose 
sight of the underlying probability space, we shall denote by P-^ (resp., E''"^ ) 
the probability measure (resp., the corresponding expectation operator) on 
this space. Finally, 9K(x) will denote the collection of all such admissible 
systems or "models" with initial configuration x = (xi,...,x„)' S (0,oo)"'. 
We shall think of 9Jl(x) as a meta-model, a collection of admissible models, 
and of the collection M = {9?^(x)}xe(o,oo)" a-s a "family of meta-models." 

The interpretation is that the components of the driving Brownian mo- 
tion W{-) represent the independent factors of the resulting model; the en- 
tries of the matrix a{t,X) are the local volatility rates of the asset capital- 
ization vector process X(-) at time t; the entries of the matrix a{t,X) as 
in (2.3) represent the local covariance rates; whereas the components of the 
vector Ti){t,X) are the market price of risk (also called relative risk) rates 
prevalent at time t. In particular, 

(2.5) j3(t, ro) := cr{t, m)'&{t, m) = y/a{t,to)'&{t, ro) 

is, in the notation of (2.3), the vector of mean rates of return for the various 
assets at time t, when the equations of (2.4) are cast in the more familiar 
form 

(2.6) dXi{t) = Xi{t)(^p^{t,X)dt + J2^aUt,X)dW,{t)y i = l,...,n. 

The integrability condition of (2.2) guarantees that the process X(-) takes 
values in (0,oo)", P-a.s.; it implies also that the exponential process 

(2.7) L(t):=exp|-^ i?'(s,X)dVF(s)-^^ ||T9(s,X)f dsj, < t < cx) 

is well defined and a strictly positive local martingale, thus also a super- 
martingale. 

This process L(-) plays the role of a state-price-density or "deflator" in 
the present context. Just as in our earlier works, Fernholz and Karatzas 
(2010a, 2010b) as well as Ruf (2011) — mostly in a Markovian context, and 
without model uncertainty — an important feature of this subject is that L(-) 
has to be allowed to be a strict local martingale, that is, that K^{L{T)) < 1 
be allowed to hold for some, if not all, T G (0,oo). 

2.2. Supermartingale problems. Constraint (1.2) brings us in the realm 
of the Krylov (1989, 2002) approach, which studies stochastic differential 
equations with constraints on the drift and diffusion coefficients in terms of 
"supermartingale problems." In particular, Theorem 2.2 of Krylov (2002) 
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shows that, under a suitable regularity condition on the family of sets IK 
in (1.1), solving stochastic equation (2.4) subject to the requirements of (2.2) 
can be cast supermartingale problem, as follows. 

Consider the nonlinear partial differential operator associated with (2.6), 
(2.5), namely 

£/(z) =F{D^f{z),Df{z),z) 

n n n \ 

ziZjQijQij + zihiPi , 

h=Jae 

(Q,p,z) gS" xM" X (0,oo)". 

The supermartingale problem is to find a probability measure P = P-^ on 
the filtered measurable space (Q,J^), F = {J^(t)}o<t<oo; under which J£(-) 
takes values in (0,00)" a.s., and the process 



u{t,X{t))- (^^{s,X{s)) + Cu{X{s))^ds, 



< t < 00 



is a local supermartingale for every u : (0, 00) x (0, 00)" — >■ M of class C^'"^ with 
compact support. 

The regularity condition on the family IC in (1.1) that we alluded to 
earlier, mandates that the function 

(2.9) F{Q,p, •) in (2.8) is Borel measurable, for every {Q,p) G S" x M". 

If, in addition, the family of sets IfC in (1.1) satisfies the linear growth con- 
dition 

(n n n \ 1/^ 

y,y2yiyj(^^j + yzivibif <c{i+ m) 
^=l ^=l i=l J 

b=^e 

and the upper-semicontinuity condition 

(2.11) limsup F{Q,p,z)<F{Q,p,y) V(Q,p) G x 

[0,oo)"9z-i>y 

for every y G &n and some real constant C > 0, then Theorem 3.2 in Krylov 
(2002) shows that the family {F^}m 6S!K(x) is convex and sequentially com- 
pact in the topology of vague convergence of probability measures. 

2.3. Markovian admissible systems. We shall also consider the subcollec- 
tion OJl^,(x) C 9JT(x) of Markovian admissible systems, for which the func- 
tionals ??(•,•) and a(-,-) as in (2.2)-(2.4) are given as 

(2.12) i^{t,X) = e{t,X{t)), a(t,X) = a(t,X(t)), 
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with measurable functions 6 : [0, oo) x (0, oo)" — ^ M", a : [0, oo) x (0, oo)" — ^ S" 
that satisfy 

(2.13) {e{t, z),a{t, z)) G /C(z) V(t, z) G [0, oo) x (0, oo)". 

Under condition (2.10) it follows then from so-called Markovian selection 
results [Krylov (1973); Stroock and Varadhan (1979), Chapter 12; Ethier 
and Kurtz (1986), Section 4.5] that the state process X{-) of (2.4) can be 
assumed to be (strongly) Markovian under P-^ , ^A G VJl^, (x) . We shall make 
this selection whenever admissible systems in 9J{*(x) are invoked. 

3. Investment rules. Consider now an investor who is "small" in the 
sense that his actions do not affect market prices. He starts with initial for- 
tune V > and uses a rule that invests a proportion Wi(t) = Ili{t,X) of cur- 
rent wealth in the ith asset of the equity market, for any given time t £ [0, oo) 
and alH = 1, . . . , n; the remaining proportion Wo{t) := 1 — "^^^i Wi{t) is held 
in cash (equivalently, in a zero-interest money market). Here n:[0, oo) x 
2IJ — > M" is a progressively measurable functional assumed to satisfy the 
requirement 



for every continuous function tt) : [0,oo) — )• (0,oo)". (Thus, the requirement 
(3.1) will be in force under all admissible systems.) We shall denote through- 
out by *p the collection of all such (nonanticipative) investment rules, and 
by the sub-collection of all Markovian investment rules, that is, those 
that can be expressed as Wi{t) = TTi{t, X(t)), < t < oo, i = 1, . . . , n for some 
measurable function vr : [0, oo) x (0, oo)*^ M". 

An investment rule is called bounded, if the functional H is bounded uni- 
formly on [0,oo) X 2IJ; for a bounded investment rule, requirement (3.1) is 
satisfied automatically, thanks to (2.2). An investment rule is called portfolio 
if the functional H satisfies Y17=i^i = 1 on [0,oo) x SB, and a portfolio is 
called long-only if Hi > 0, ... , n.„ > also hold on this domain. A long-only 
portfolio is clearly bounded. 

Given an initial wealth v G (0, oo), an investment rule H G ^ and an admis- 
sible model M G 9Jt(x), the resulting wealth process Z{- ) = Z^'n(-) satisfies 
the dynamics 



(3.2) ^ = ^m,X) = U'{t, X)a{t, X) mt, X) dt + dW{t)] 



and the initial condition Z{0) = v. In conjunction with (2.7) in the differen- 
tial form 



(3.1) 




(|n'(t, m)a{t, Ki)^(t, m)\ + n'(t, to)a{t, m)U(t, ra)) di < oo 



VTg (0,oo) 



n 



1=1 



(3.3) 



dL{t) = -L{t){'d{t,X)ydW{t) 
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and the product rule of tfie stocliastic calculus, this gives 

L{t)Z'''^{t) =v+ f L{s)Z'''^{s){a'{s, X)n(s, X) - X))' dW{s), 
(3.4) 

< t < oo. 

For any initial configuration x = G (0,oo)", initial wealth v G 

(0, oo), investment rule 11 S *P and admissible model M. € 9Jl(x), the product 
L{-)Z'"'^{-) is therefore under P-^ a continuous, positive local martingale, 
thus also a supermartingale. Once again, it is important that this process 
be allowed to be a strict local martingale. 

3.1. The market portfolio. The choice of Markovian investment rule m G 
given by 

mi{t, z) = — ' — , z = 1, . . . , n, t G [0, oo), z G (0, oo)" 

^1 + ■ ■ ■ + Zn 

leads to the long-only market portfolio 

^i{t) = X{t)/X{t) 

(3.5) 

with weights ^i{t) = Xi{t)/X{t), i = I, . . . ,n, 0<t<oo 

in the notation of (2.1). It follows from (3.2) that investing according to this 
portfolio amounts to owning the entire market, in proportion to the initial 
wealth, Z^'''^{-) = vX{-)/X{0). 

3.2. Ramifications. Reading (3.4) for the market portfolio of (3.5), and 
recalling (2.1), leads to 

L{t)X{t) = X{0)+ [ L{s){a'{s,X)X{s) - ^s,X)X{s)y dW{s), 
(3.6) ^° 

< t < OO 

or equivalents d{L{t)X{t)) = -L{t)X{t)(&{t,X)y dW{t), where 
TT/ N n/ N a'(t,xv)tv(t) 

^ ' ^ ^ tt)i(t) + --- + tt)n(t) 

(3-7) ^ 



satisfies / ||??(t, tt))|p dt < oo 
Jo 



for all (t, ro) G [0,oo) x 22J, thanks to (2.2). With this notation, it follows 
from (3.6) that 

L{.)X{-) = {xi + --- + Xn)-expl- f {d{t,X)ydW{t)-^ f Mt,X)fdt\. 
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On the strength of the integrabiUty condition in (3.7), the Dambis-Dubins- 
Schwartz representation [e.g., Karatzas and Shreve (1991), page 174] of the 
* -local martingale 



N{.):= / {^t,X)ydW{t) 
Jo 

with quadratic variation {N){-) = / ||T?(t, jC)|p dt < oo 

Jo 

gives 

(3.8) L(r)X(r) = (xi + --- + x„)-e^(")-("/2)|^^^^^^^^, 0<r<oo, 

where B{-) is one-dimensional, standard Brownian motion under P-^. 
Whereas the equations of (2.4) can be written as 

(3.9) dX,{t) = Xi{t) f dt + aUt, X) dW.(t) J 
for i = 1, . . . ,n, with 

(3.10) W{-):=W{-)+ [ d{t,X)dt. 

Jo 

We then have the representation 

A(-) := -^^^ = expM (^?(t,X))MT^(t) + i^ |Kt,X) 
(3.11) 

= exp I ^ {d{t, X))' dW{t) - i ^ mt, X) f di I 

for the normalized reciprocal of the deflated total market capitalization, and 

(3.12) dfii{t) = fXi{t){ci-fi{t))'a{t,X)dW{t), i = l,...,n 

for the dynamics market weights in (3.5); here ti is the ith unit vector in M". 

4. Optimal arbitrage relative to the market. Let us consider now the 
smallest proportion 

u(T, x) = inf {r > : 3n^ G ^, s.t. P"^ (Z'-^W.nr (T) > X(r)) = 1, 

(4.1) 

VA^ G 9Jt(x)} 

of the initial total market capitalization X (0) = xi + ■ ■ ■ + Xn which allows 
the small investor, starting with initial capital u(T, x)X(O) and through 
judicious choice of investment rule in the class to match or exceed the 
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performance of the market portfolio over the time- horizon [0,T], and to 
do this with P-^-probabihty one, under any model A4 G 9Jt(x) that might 
materialize. We shall refer to u(-, •) of (4.1) as the arbitrage function for the 
family of meta-models M = {9?t(x)}xg(o,oo)" ; think of it as a version of 
the arbitrage function studied in Fernholz and Karatzas (2010a) which is 
"robust" with respect to M. 

The quantity of (4.1) is strictly positive; see Proposition 1 below and 
the discussion following it. On the other hand, the set of (4.1) contains the 
number r = 1 , so clearly 

0<u(r,x) < 1. 

If u(r, x) < 1, then for every r E (u(r, x),l) — and even for r = u(T, x) 
when the infimum in (4.1) is attained, as indeed it is in the context of 
Theorem 1 below — there exists an investment rule 11^ S *p such that 

^x(o),nw^N > }ix(T) > X(T) = Z^(°)''^(r), P-^-a.s. 
r 

holds for every M € 9JT(x). In other words, the investment rule 11^ leads 
then to strong arbitrage relative to the market portfolio in the terminology 
of Fernholz and Karatzas (2009) — here with the extra feature that such 
arbitrage is now robust, that is, holds under any possible admissible system or 
"model" that might materialize. If, on the other hand, u(r, x) = 1, then such 
outperformance of (equivalently, strong arbitrage relative to) the market is 
just not possible over all meta-models M G 9Jt(x). In either case, the highest 
return on investment relative to the market 

b(r, x) := sup{6 > : 3n G *p, s.t. p^(z^(o),n(^) > bX{T)) = 1, 

yM G 93T(x)}, 

achievable using (nonanticipative) investment rules, is given as b(T, x) = 
l/u(T,x)>l. 

Remark 1. Instances of u(T, x) < 1 occur, when there exists a constant 
C > such that either 

(n n n \ 

— zm^ yy- — >^ 
1 = 1 1 = 1 J = l / 

or 

(4.3) ('(fiiiifiOi^y ij^f (yau--yya,]>c 

holds for every z G (0,c«)". See the survey paper Fernholz and Karatzas 
(2009), Examples 11.1 and 11.2 [as well as Fernholz and Karatzas (2005), 
Fernholz, Karatzas and Kardaras (2005) for additional examples]. 
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Proposition 1. The quantity of (4-1) satisfies 

u(T,x) >^>(r,x) >0 

(4.4) 

where «>(r,x) := sup ' L v ; v ;j \ 



M&<SX{x) \ H \- Xn J 

Furthermore, under conditions (2.9)-(2.11), there exists an admissible sys- 
tem Mo G 97t(x) such that 

Proof. Take an arbitrary element r > of the set on the right-hand side 
of (4.1) and an arbitrary admissible system M € 9Jt(x). There exists then 
an investment rule 11^ G «p with the inequality Z''^(°)'"'-(r) > X{T) vahd 
P-^-a.s. On the strength of (2.7) and (3.4), the process L{^Z''^^^^'^^ {■) is 
a P-^-supermartingale; thus (3.8) and (3.7) lead to 

r{xi + ---+Xn)= rX(0) > E^-^ [L{T)Z'-^'^^'^^^- (T)] 

(4.6) 

>E^-^[L(r)x(r)] >o. 

The inequality u(r, x) > $(r, x) in (4.4) follows now from the arbitrariness 
of r > and M G 9JI(x). The existence of an admissible system Mo G S[R(x) 
that satisfies (4.5) follows from Theorem 3.4 in Krylov (2002), in conjunction 
with the dynamics of (2.4) and (3.3). □ 

Although strong arbitrage relative to the market may exist within the 
framework of the models M G OJt(x) studied here (cf. Remark 1), the exis- 
tence of a strictly positive supermartingale deflator process L(-) as in (2.7) 
proscribes scalable arbitrage opportunities, also known as Unbounded Profits 
with Bounded Risk (UPBR); this is reflected in the inequality u(T, x) > 
of (4.4). We refer the reader to Delbaen and Schachermayer (1995b) for 
the origin of the resulting NUPBR concept, and to Karatzas and Kardaras 
(2007) for an elaboration of this point in a different context, namely, the 
existence and properties of the numeraire portfolio. 

Finally, let us write (4.4) as 

(4.7) «>(r,x)= sup UA,(r,x) where UA,(r,x):= ^^ [L{T)X{T)\ ^ 
xean(x) ^i^ rXn 

We have for this quantity the interpretation 

UM {T, x) = inf {r > : G «P, s.t. (Z-^^W-n. (5-) > x{T)) = 1} 
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as the "arbitrage function for the model M € 9K(x)," at least when the 
matrix a{t, w) in invertible for every (t , lu) G (0, oo) X 2B and when (P^ , F)- 
martingales can be represented as stochastic integrals with respect to the 
Brownian motion W{-) in (2.4). 

5. A fully nonlinear PDI. Consider now a continuous function U : [0, oo) x 
(0,oo)"^ (0,oo) with 

(5.1) C/(0,z) = l, zG(0,oo)", 

which is of class C^'^ on (0,oo) x (0,00)" and satisfies on this domain the 
fully nonlinear partial differential inequality (PDI) 

1=1 7=1 ^ i I I n/ 

(5.2) 

Vo e A{z). 

We shall denote by U the collection of all such continuous functions U : [0, 00) x 
(0, 00)" — ;> (0, 00) which are of class C^''^ on (0, 00) x (0, oo)'^ and satisfy (5.1) 
and (5.2). This collection U is nonempty, since we can take U{-,-) = 1; how- 
ever, l/( need not contain only one element. 

Let us fix an initial configuration x G (0,00)" and consider any admissible 
system Ai E lUl(x). Applying Ito's rule to the process 

(5.3) E{t):=X{t)L{t)U{T -t,X{t)), 0<t<T 

in conjunction with (3.6) and (2.4), we obtain its P-^-semimartingale de- 
composition as 



+ ^[i?^(t, X) - U{T - t, X{t))Mt, ^)] dW^it). 
Here we have used the notation of (3.7), and have set 

n 

(5.5) R,(t,T) ■.= Y,Mt)DiU{T -t,X{t))ai,{t,X), 



1=1 

n n 



A(i, ^):=\Y.iZ Mt)Xjit)a,j (t, X)Df^U{T - t, X(t)) 



2 

i=i j=i 



i — 1 \j — 1 / 



dU 

— (T-t,X(t)). 
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From the inequality of (5.2), coupled with the fact that a{t,X) € A{3C{t)) 
holds for all < t < oo, this last expression is clearly not positive. As a result, 
the positive process H(-) of (5.3) is a P-^-supermartingale, namely, 

L{t)X(t)U{T - t, X(t)) = E{t) > E^^ [E{T)\T{t)] 

(5.7) 



E^'^[L{T)X{T)\J^{t)] 



holds P-^-a.s., yM G 9Jt(x) and < t < T; in particular, 

(xi + • • • + x„)C/(r, x) = S(0) > E^^ [E{T)] 

= E^^[L{T)X{T)] VA^G5Jt(x). 

With the notation of (4.4), we obtain in this manner the following analog of 
the inequality in Proposition 1: 

(5.8) [/(r,x)>a>(T,x). 

Digging in this same spot, just a bit deeper, leads to our next result; this is 
very much in the spirit of Theorem 5 in Fleming and Vermes (1989) and of 
Section II.2 in Lions (1984). 

Proposition 2. For every horizon T S (0, oo), initial configuration x E 
(0,oo)" and function U : [0,oo) x (0,oo)" — )• (0,oo) in the collection U, we 
have the inequality 

(5.9) c/(r,x)>u(r,x)>$(r,x). 

Furthermore, the Markovian investment rule G generated by this func- 
tion U through 

^f(t,z) :=ZiD,logU{T-t,z) 

(5.10) 

+ , (t,z)G[0,r]x(0,oo)- 

Zi -r • • • -\- Z^i 

for each i = 1, . . . ,n, satisfies for every admissible system Ai G 9Jt(x) the 
inequality 

(5.11) Z^('^'-)^(0)'-^(r)>X(r), ¥^-a.s. 

Proof. For a fixed initial configuration x G (0, oo)", an arbitrary admis- 
sible model G 9Jt(x) and any function U £U, let us recall the notation 
of (5.3) and re-cast the dynamics of (5.4) as 

(5.12) dE{t) = E{t) (j2^ ^,{t, X) dW,{t) - dC{t)j . 
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Here by virtue of (5.4), (5.5) and (3.7) we have written 

^,{t, X) := J2 ^-(^' ^) A log UiT - t, X{t)) + ^) 

(5.13) 

for ly = 1, . . . ,n and have introduced in the notation of (5.6) the continuous, 
increasing process 

(5.14) C{t):= / \i~^^'''S\, ds, 0<t<r. 



JO U{T-s,X{s)) 

The expression of (5.13) suggests considering the Markovian investment 
rule vr^ E as in (5.10); then we cast the expression of (5.13) as 

n 

^,{t, X) = ^ aUt, X)7r^(t, X{t)) - Mt, 
1=1 

On the strength of (3.4), the value process generated by this investment 
rule TT^ starting with initial wealth ^ := [/(T, x)X(O) = H(0), satisfies the 
equation 

n 

d(L(t)Z«'-'' (t)) = {L{t)Z<^^^" (t)) "fuit, X) dW,{t). 

v=l 

Juxtaposing this to (5.12), and using the positivity of H(-) along with the 
nonnegativity and nondecrease of C(-), we obtain the P^-a.s. comparison 
L(-)Z«'-''(-) >H(-), thus 

(5.15) Z^^^" {t)>X{t)U{T -t,X{t)), 0<t<T. 

With t = T this leads to (5.11), in conjunction with (5.1). We conclude 
from (5.11) that the number [/(T, x) > belongs to the set on the right- 
hand side of (4.1), and the first comparison in (5.9) follows; the second is 
just a restatement of (4.4). □ 

Corollary. Suppose that the function <!>(•,•) of (4-4) belongs to the 
collection U. Then <&(•,•) is the smallest element ofU; the infimum in (4-1) 
is attained; we can take U = ^ in (5.11), (5.10); and the inequality in (5.9) 
holds as equality, that is, <&(•,•) coincides with the arbitrage function 

u(r,x) = $(r,x) v(r,x) G (o,oo) x (o,oo)". 

Interpretation: Imagine that the smah investor is a manager who invests 
for a pension fund and tries to track or exceed the performance of an index 
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(the market portfolio) over a finite time-horizon. He has to do this in the face 
of uncertainty about the characteristics of the market, including its covari- 
ance and price-of-risk structure, so he acts with extreme prudence and tries 
to protect his chents against the most adverse market configurations imag- 
inable [the range of such configurations is captured by the constraints (1.2), 
(1.1)]. If such adverse circumstances do not materiahze his strategy gener- 
ates a surplus, captured here by the increasing process C(-) of (5.14) with 
U = $ = u, which can then be returned to the (participants in the) fund. 
We are borrowing and adapting this interpretation from Lyons (1995). 

Similarly, the Markovian investment rule tt^ G generated by the func- 
tion [/ = <I> = u in (5.10), (1.3) implements the best possible outperformance 
of the market portfolio, as in (5.11). 

6. A verification result. For the purposes of this section we shall impose 
the following growth condition on the family A = {^(y)}ye6„ of subsets 
of in (1.6), (1.1): there exists a constant C G (0,oo), such that for all 
y G 6n we have 

(6.1) sup (max , ^f^'^^f' . ) < C{1 + ||y||). 

aeA{y) yi^^^j^" {yi + --- + yn)J 

We shall also need the following strong ellipticity condition, which mandates 
that for every nonempty, compact subset K of (0,oo)", there exists a real 
constant A = Ak > such that 

(6.2) inf( inf jX^X^e^e.a., ) ) >AK||^f G M" 

Assumption A. There exist a continuous function a : (0, oo)" — )• S*^, 
a C^-function : (0, oo)"' — > M, and a continuous square root s(-) of a(-), 
namely a(-) = s(-)s'(-) such that, with the vector-valued function 6{-) = 
(0i(-),...,0n(-)y defined by 

n 

(6.3) 6^{z) ■.= '^ZjSj„{z)DjH{z), iy = l,...,n, 

i=i 

condition (2.13) is satisfied, whereas the system of stochastic differential 
equations 

n 

dXi(t) = ^Si„(S(t))[(«F„(t) + e„(3C(tj)dt], 

(0.4) 

Xi[0) = Xi > 0,i = 1, . . . ,n 
has a solution in which the state process X(-) takes values in (0,oo)". 
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A bit more precisely, this assumption posits the existence of a Markovian 
admissible system ^Ao £ S!Jt*(x) consisting of a filtered probability space 
(r2,J^,P), F = {J^(t)}o<t<oo and of two continuous, adapted process j£(-) 
and W{-) on it, such that under the probability measure P = P-^° the pro- 
cess W{-) is n-dimensional Brownian Motion, the process X(-) takes values 
in (0,oo)"' a.s. and (2.4) holds with 'd,y{t,X) = 6^{X{t)) as in (6.3), and with 
a"jjy(t, X) = Si^{X{t)), < t < oo (1 < i, < n). The system of equations (6.4) 
can be cast equivalently as 



(6.5) 



^Si,(X(t))dVF,(t) 



+ ^aij{X{t))Xj{t)DjH{Xit))j dt 



Assumption B. In the notation of the previous paragraph and under 
the condition 

n n 

(6-6) 5^5^Zi|si,(z)||0,(z)| <C(l + ||z||) VzG(0,oo)", 

i=l v=l 

we define on (0, oo)" the continuous functions g{z) := e~^^^^ Y17=i ^^"^ 
A;(z) := {l/2)i2'l=i^Z'j=i^3{^)[D'fjH{x) + DiH{2)DjH{'z)\ and assume that 
the function 



G(r,x) :=E 



c/(X(T))exp<^ / k{X{t))dt 



(r,x)G[0,oo)" 



is continuous on [0,oo) x (0, oo) and of class C ' on (0, oo) x (0,oo). 

Sufficient conditions for Assumptions A, B to hold are provided in Fern- 
holz and Karatzas (2010a), Sections 8 and 9. It is also shown there, that we 
have the P-'^°-martingale property 

(6.7) wF^''[X{T)L{T)\F{t)]=X{t)L{t)-T{T -t,X{t)), 0<t<T 
for the function 

(6.8) r(T,z):=G(r,z)/r7(z), (r, z) G [0, oo) x (0, oo)". 

This function is of class C^'^ on (0,oo) x (0,oo) and satisfies the initial con- 
dition r(0,-) = 1 on (0,00)" as well as the linear second-order parabolic 
equation 

Ar(T,z) 



(6.9) 



^(^, z) = ^ ^ ZiZ,^i,{i) [ -D%Y[j, z) + 

i=\ 1=1 



Z\^ V Zn, 

(t,z) G (0,00) X (0,00)". 
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Proposition 3 (Verification argument). Under the Assumptions A, B 
and the conditions (6.1), (6.2), suppose that the functions a(-) and r(T,-) 
satisfy the inequality 

i=l j=l ^ i I I n/ 

(6.10) 

>±± za,a., (ll)^.r(r, z) + Va G ^(z) 

for every {t,z) G (0,cx3) x (0,oo)". Then, in the notation of (4.1)-(4-'^), we 
have: 

u(T, x) = ^T, x) = r(r, x) = UM, (T, x), V(r, x) G (0, oo) X (0, oo)" 

for the Markovian admissible system M = Mo G S!?t(x) posited in Assump- 
tion A; the conclusions of Proposition 2 and its corollary for U = as well 
as the P-^-a.s. comparison 

L{t)X{t)-u{T-t,X{t))>E^^[L{T)X{T)\T{t)], 0<t<T, 
which holds for every Ai G 9H(x) and as equality for A4 = A^o G 9?t(x). 

Proof. Under condition (6.10) the function r(-,-) belongs to the col- 
lection U, as (5.2) is satisfied with ?7 = F on the strength of (6.9) and (6.10); 
thus, we deduce r(r, x) > $(T, x) from (5.8). On the other hand, equal- 
ity (6.7) with t = 0, and the definition of <^>(T, x) in (4.4), give 

Xi H \-Xn 

so the equality r(r, x) = $(T, x) follows. In other words, we identify A^o as 
a Markovian admissible system that satisfies (4.5) and attains the supremum 
in (4.4). The remaining claims come from Proposition 2 and its Corollary 
[in particular, from reading (5.9) with U = T] and from (5.7), (6.7). □ 

Remark 2. Proposition 3 holds under conditions weaker than those 
imposed in Assumptions A and B above, at the "expense" of a certain lo- 
calization. More precisely, one posits the existence of locally bounded and 
locally Lipschitz functions Sjjy(-) and Oi^{-) (1 < < ra) for which (2.13), 
(6.1), (6.2) and (6.6) are satisfied with a(-) =s(-)s'(-), and for which there 
exists a Markovian admissible system Mo G 9Jt*(x) whose state process X(-) 
in (6.4) is, under F"^", a strong Markov process with values in (0,oo)" a.s. 
Using results from the theory of stochastic flows [Kunita (1990), Protter 
(2004)] and from parabolic partial differential equations [Janson and Tysk 
(2006), Ekstrom and Tysk (2009)], Theorem 2 in Ruf (2011) shows that the 
function r(-,-) is then of class C^'^ locally on (0,oo) x (0, oo)", and solves 
there equation (6.9). 
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7. Maximizing the probability of containment. We have now gone as far 
as we could without having to specify the nature of our filtered measurable 
space ¥ = {-7^(t)}o<t<oo- To proceed further, we shall need to choose 

this space carefully. 

We shall take as our sample space the set ft of right-continuous paths 
a; : [0, oo) — S„ U {A} . Here A is an additional "absorbing point" ; paths stay 
at A once they get there, that is, after T(w) = inf{t > 0|a;(i) = A}, and are 
continuous on (0, ^(a;)); we are employing here, and throughout this work, 
the usual convention inf = oo. We also select /C(A) = {(O,0„xn)} and 
^(A) = {0 

nxn}) where O^xn is the zero matrix. With P{t) := a{ui{s),0< 
s<t), the filtration = {P{t)}o<t< oo is a standard system in the terminol- 
ogy of Parthasarathy (1967). This means that each (^},J-^{t)) is isomorphic 
to the Borel cr-algebra on some Polish space, and that for any decreasing 
sequence {Aj}j^fq where each Aj is an atom of the corresponding J^{tj), 
for some increasing sequence {ijjjeN C [0,oo), we have fljeN^i ^ i^^^ 
Appendix in Follmer (1972), as well as Meyer (1972) and Follmer (1973)]. 

With all this in place we take (O, J^),F = {-7^(t)}o<t<oo as our filtered 
measurable space, where 



An admissible system M G 5Jt(x), x G (0,oo)" defined as in Section 2 con- 
sists of this filtered measurable space (O, J^),F = {J^{t)}o<t<oo, of a proba- 
bility measure P-^ on it, of an n-dimensional Brownian motion W{-) on the 
resulting probability space and of the coordinate mapping process X{t,ui) = 
u){t),0 <t<oo which is assumed to satisfy (2.4), (2.2) and to take values in 
(0,oo)", P-a.s. We shah take 

r := inf{t > 0|A(t) = 0} = mi{t > 0\L{t)X{t) = oo} 

in the notation of (3.11), (2.7) and (2.1), and note P-^(r < oo) = 0. 

7.1. The Follmer exit measure. With this setup, there exists a probabil- 
ity measure Q on (i7, J^), such that 

(7.1) dP-^ = A(r)dQ holds on each J^(r), re(0,oo); 

we express this property (7.1) by writing P-'^ ^ Q. Under the measure Q, the 
process W{-) of (3.10) is Brownian motion; whereas the processes /ii(-), . . . , 
/!„(•) and A(-) of (3.12), (3.11) in Section 3.1 are nonnegative Q-martingales. 

The "absorbing state" A acts here as a proxy for P-'^-nuU sets to which the 
new measure Q may assign positive mass; the possible existence of such sets 
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makes it important that the filtration F be "pure," that is, not completed 
by P-^-null sets. This probability measure Q satisfies 

^"■^*^'^-<^" -K- Id/Am) ^ 



xi-\ h x„ 

■ > T) VT G [0, oo) 



(7.2) 
and 

(7.3) r = inf|i>0 ||^(s,X)f ds = oo|, Q-a.s. 

We also have Q-a.e. on {T < T < oo} 

L{T + u)X{T + u) = oo \/u>0 

and 

\\d{t,X)fdt< r \\d{t,X)fdt = oo. 



I 

Jo 



/o Jo 
Whereas, Q-a.e. on {T = oo}, we have 

L(r)X(r)<oo, / ||T?(t,X)f dt <oo; VrG[0,oo). 
Jo 

We deduce from (7.2) that the arbitrage function of (4.7) for the model 
A4 £ 9?I(x) is given by the probability of ^^containment" under the mea- 
sure Q, namely, the Q-probability that the process X(-), started at x S 
(0,oo)", stays in (0,oo)" throughout the time-horizon [0,r]. 

At this point we shall impose the following requirements on K = {^(y)}yes„ i 
the family of compact, convex subsets of M*^ x S" in (1.1), (1.2): there exists 
a constant < C < oo, such that for all y G ©.„ we have the strengthening 

(n n \ 
"^"^yiVjaij <C{yi-\ \-ynf 
i=i i=i J 

of the growth condition in (6.1), as well as the "shear" condition 
(7.5) sup 

(e,a)G/C{y) 



+ 



< C. 



l-hTr(a) J Vl + II^P 
Then the following identity holds Q-a.s.: 

(7.6) r= min T; where 71 := inf{t > 0|Xi(t) = 0}. 

l<i<n 

For justification of the claims made in this subsection, we refer to Section 7 
in Fernholz and Karatzas (2010a), as well as Delbaen and Schachermayer 
(1995a), Pal and Protter (2010) and Ruf (2011)— in addition, of course, to 
the seminal work by Follmer (1972, 1973). 
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The special structure of the filtered measurable space {Q,T),¥ = 
{J^(i)}o<t<oo that we selected in this section is indispensable for this con- 
struction and for the representation (7.2); whereas the inequality < 
C(l + Tr(a)), V(^, a) E K,{y),y E S.„ from condition (7.5) is important for 
establishing the representation of (7.6). 

7.2. Auxiliary admissible systems. Let us fix then an initial configuration 
x= (xi, . . . jXn)' E (0,oo)" and denote by 5'I(x) the collection of stochas- 
tic systems M that consist of the filtered measurable space (il, J^),F = 
{J^{t)}o<t<oo, of a probability measure Q = Q^, of an M"-valued Brown- 
ian motion W(-) under Q, and of the coordinate mapping process X{t,ui) = 
uj{t),{t,uj) E [0,oo) X which satisfies Q-a.s. the system of the stochastic 
equations in (3.9) 

dx.(t) = dt + E aut, X) dwAt)j 

(7.7) ^ x.it) t .ut. X) ( dm) + 1 ^) , 

i = 1, . . . , n. 

Here the elements aii, : [0, oo) x — M, 1 < i,u < n oi the matrix cj(-, •) = 
{(^iui'^ ')}i<i,u<n are progressively measurable functionals that satisfy, in the 
notation of (1.6), 

(7.8) (T{t,uj)(T'{t,uj)=:a{t,oj) (^A{uj{t)) V(t,a;) E [0, oo) x 

As in Section 2.3, we shall denote by 9T^,(x) the subcollection of ^(x) that 
consists of Markovian auxiliary admissible systems, namely, those for which 
the equations of (7.7) are satisfied with a{t,X) =a(t,X(t)) and a(t,X) = 
s{t,X{t)), < t < oo and with measurable functions a: [0,oo) x (3„ — t- 
and s:[0,oo) x 6„ — )-£(M";M'^) that satisfy the condition s(i,y)s'(t,y) = 
a(t,y) E A{y), V(t,y) E [0,oo) x We invoke the same Markovian selec- 
tion results as in Section 2.3, to ensure that the process X(-) is strongly 
Markovian under any given Q^, Af E 9T*(x). 

By analogy with (7.6), we consider 

(7.9) r(a;) := min 7;(a;) with 7I(a;) = inf{i > 0|wi(t) = 0}. 

l<i<n 

Then for every oj & {T < oo} we have 

(7.10) / Tr{a{t,oj))dt< / Ti {a{t, lo)) dt = oo \fO <T <T{uj); 
Jo Jo 

whereas Tr(a(t, oo)) dt < oo, < T < oo holds for every w E {T = oo}. 
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Remark 3. As in Section 2.2, solving the stochastic equation (7.7) sub- 
ject to condition (7.8) amounts to requiring that the process 

u{t, X{t)) - (^(s, ^{s)) + Cu{X{s))^ ds, < t < oo 

be a local supermartingale, for every continuous u: (0,oo) x ©„ — t- M which 
is of class C^'^ on (0,oo) x (OjOo)" and has compact support; here C is the 
nonlinear second-order partial differential operator in (1.5). 

Remark 4. The total capitalization process X{-) = Xi{-) + • • • + Xn{-) 
satisfies, by virtue of (7.7), the equation 

dX{t) =X{t)[dN{t) + d{N){t)l 



v=\-'^ \l=\ / 

Under the measure Q, the process A^(-) is a continuous local martingale with 
quadratic variation 

n n „i 

(iV)(t) = VV / Xi[s)ai,[s,X)X,[s){X^{s)^---^Xn[s)r''ds<Ct 

from (7.4), so the total capitalization process 
X{t) = X{0) ■ e^W+(V2)W(-) = ^(0) . e5M+K2)|^^^^^^^^^ < t < oo 

takes values in (0,oo), Q-a.e.; here -B(-) is a one-dimensional Q-Brownian 
motion. This is in accordance with our selection of the punctured nonnega- 
tive orthant &n in (1.1) as the state-space for the process X(-) under Q. 

Under this measure, the relative weights /ii(-) = Xi{-)/X{-), i = 1, . . . ,n 
are nonnegative local martingales and supermartingales, in accordance with 
(3.12), and since X^iLi/^j(") = 1 these processes are bounded, so they are 
actually martingales. Once any one of the processes ^i(-)i • • • > Xn{-) [i.e., any 
one of the processes /xi(-), //„(•)] becomes zero, it stays at zero forever; 
of course, not all of them can vanish at the same time. 

In Section 7.1 we started with an arbitrary admissible system A4 G 9Jt(x) 
and produced an "auxiliary" admissible system M G ^(x), for which the 
property (7.2) holds. Thus, for every (T, x) G (0,oo) x (0,oo)" we deduce 

g(r,x):= sup Q^iT>T) 

(7.11) 

'E^^ [L{T)X{T)]\ 

M&M{x)\ H \-Xr. 



> sup ( ^ , ; l=<I>(T,x). 
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7.3. Preparatory steps. We suppose from now onwards that, for every 
progressively measurable functional a : [0, oo) x — t- S" which satisfies 

(7.12) a{t,oj)eA{uj{t)) for all (t,w) G [0,oo) X 1^, 

we can select a progressively measurable functional : [0, oo) x 17 — )■ M" with 

(7.13) {^{t,uj),a{t,u})) elC{uj{t)) y{t,uj) e[0,oo) X n 

[see the "measurable selection" results in Chapter 7 of Bertsekas and Shreve 
(1978)]. We introduce now the functional 

(7.14) d{t,u}) := i3{t,uj) - a'{t,uj)uj{t)/{uji{t) + • • • + a;„(t)) 
as in (3.7) and also, by analogy with (7.3), the stopping rule 

(7.15) r(w) :=inf|t >0 ^ ||??(s, w)|P ds = oo| 

[cf. Levental and Skorohod (1995), where stopping rules of this type also 
play very important roles in the study of arbitrage]. 

We recall now (7.10); on the strength of the requirement (7.5), this gives 

(7.16) / \\^{t,uj)fdt< / \\^{t,u})fdt = oo, 0<r<r(a;) 



JO JO 

for every oj G {T < oo}, and Jq \\-d{t,uj)\\'^ dt < oo, VT G [0, oo) for every 
u G {T = oo}. In conjunction with (7.4), we obtain from (7.16) that 

T ^ rf{uj) _ 

||??(t,w)f dt< / ||i9(t,w)f dt = oo, o<r<r(L^) 

JO 

holds for every uj G {T < oo}, and that \\-d{t,oj)\\'^ dt < oo, VT G [0,oo) 

holds for every a; G {T = oo}. 

We deduce for the stopping rules of (7.15) and (7.9) the identification 
f(a;)=r(a;). 

7.4. The same thread, in reverse. Let us fix now a stochastic system 
M G 9T(x) as in Section 7.2, pick a progressively measurable functional 
a: [0,oo) X Q with a{t,u) G A{uj{t)) for all (t,w) G [0,oo) x Q. and 

select a progressively measurable functional -d : [0, oo) x J] — )• as in (7.13). 
For this '&{■, •) and this M G 01(x), we define i?(-, •) by (7.14) as well as 

A(t) = expj j'(d{s, X))' dW{.s) -^1' Ms, X) IP ds 
(7.17) ° ° 

for < t < r 
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as in (3.11), and set 

(7.18) A(r + M) = for u>0 on {r<oo} 

in the notation of (7.15). The resulting process A(-) is a local martingale 
and a supermartingale under Q, and we have T{3i) = inf{t > 0|A(t) = 0}, 
Q-a.e. 

We introduce also the sequence of F-stopping rules 

S'„(a;) :=inf|t>0 ||i9(s,tj)f ds > 21ogn|, n € N, 

which satisfy lim„_j.oo t Snioj) = T{oj) and exp{i /^f"^"^^ \\'d{t,oj)f dt] < n 
(Vn G N), for every lo £ Q. From Novikov's theorem [e.g., Karatzas and 
Shreve (1991), page 198], A(- A Sn) is a uniformly integrable Q-martingale; 
in particular, E''2(A(S'„)) = 1 holds for every n G N. Thus, the recipe 

P„(yl):=E^[A(5„)-lA], AeT{Sn) 

defines a consistent sequence, or "tower," of probability measures {P^jneN 
on (r2,J^). Appealing to the results in Parthasarathy (1967), pages 140- 
143 [see also the Appendix of Follmer (1972)], we deduce the existence of 
a probability measure P on {i^,J-) such that 

(7.19) ¥{A) = Fn{A) = E^[A{Sn) • Ia] holds for every A G J^(5„),n G N. 

(Here again, the special structure imposed in this section on the filtered 
measurable space {Q,,T),¥ = {T{t)}o<t<oo is indispensable.) Therefore, for 
every T G (0, oo) we have 

P(5„ > T) = E^[A{Sn) ■ 1{5„>T}] = E^[A(r) • 1{5„>T}] 

by optional sampling, whereas monotone convergence leads to 

(7.20) P(r>r) = E^[A(T)l|r>T}]- 

The following result echoes similar themes in Cheridito, Filipovic and Yor 
(2005). 

Lemma 1. The process A(-) of (7.17), (7.18) is a Q-martingale, if and 
only if we have 

p(r<oo) = o 

[i.e., if and only if the process X(-) never hits the boundary of the orthant 
(0,oo)", F-a.s.]. 

Proof. If P(T < co) = holds, the nonnegativity of A(-) and (7.20) 
give 

1 = P(r > T) = E^[A(T)l{r>T}] < E*^[A(r)] VT G (0, oo). 
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But A(-) is a Q-supermartingale, so the reverse inequality E'^[A(T)] < A(0) = 
1 also holds. We conclude that E^[A(r)] = 1 holds for all T G (0,oo), so A(-) 
is a Q-martingale. 

If, on the other hand, A(-) is a Q-martingale, then EQ[A(r)] = 1 and (7.20) 
give 

P(r < T) = E^(A(r)) - E^{A{T)l^r>T}) = E^(A(r)l|^<T|) 
= E^(A(r)l{r<T})=0 

for every T £ [0,oo), from Optional Sampling and the fact that A(7') = 
holds Q-a.e. on {T < oo}. We conclude P(T < oo) = 0; in conjunction with 
the identification T = T and (7.9), this means that the coordinate mapping 
process never reaches the boundary of (i.e., takes values in) the strictly 
positive orthant (0,cxd)", P-a.e. □ 

When the conditions of Lemma 1 prevail, the process 

W{-) = W{-)- [ d{t,X)dt 
Jo 

is Brownian motion under the probability measure P = P''^ introduced 
in (7.19). This measure satisfies the equations of (7.2), whereas the pro- 
cess X(-) solves P-^-a.s. the system 

n 

dXi{t)=X,{t)^ai^{t,X)[dW^{t) + Mt,^)dt], X,{0) = Xi > 
i/=i 

for i = 1, . . . ,n, as in (2.4). It is then not hard to check that L(-) defined 

by (2.7) satisfies P-^-a.s. the identity L{-)X{-) = (xi H h Xn)/A{-) in 

accordance with (3.11). 

We formalize these considerations as follows. 

Assumption C. Suppose that the collection of sets IC in (1.1) satis- 
fies (7.4), (7.5) and that for any given progressively measurable functional 
a: [0,oo) X — which satisfies (7.12) and 

(7.21) / Ti'(a(t,tt)))dt <oo for all (T, ro) E [0, oo) x 2U, 
Jo 

there exists a progressively measurable functional -d : [0, oo) x i7 — )■ M" that 
satisfies the condition (7.13), thus also by virtue of (7.5) 

(7.22) [ \\^{t,n)f dt < oo for all (T,rti) G [0,oo) X 2n. 
Jo 
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The analysis of this subsection shows that, under Assumption C and start- 
ing with any initial configuration x = (xi , . . . , G (0, 00)*^ and with an ar- 
bitrary "auxiliary" admissible stochastic system AA= ((fi, J^), F,Q,X(-), 
in 9T(x) as in Section 7.2, the process A(-) of (7.18) is a Q-martingale, and 
we can construct a "primal" admissible system A4 G 9Jt(x) as in Section 2.1 
[i.e., with the canonical process X(-) taking values in (0,oo)", P-^-a.s.], for 
which (7.2) holds, and we have <C Q as in (7.1). We deduce 

^ E^^ [L{T)X{T )]\ 
ATe'oiCx) ' ' MeW-ix) V H \- Xr, 



Q{T,^)= sup (^^T>T)< sup ( ^ , , : 1 =^(T,x). 



The reverse inequality Q{T,x.) > <1>(T, x) was established in (7.11). This way, 
for every function U : [0,oo) x (0,oo)" — t- (0, 00) in the collection U, we can 
strengthen (5.9) to 

c/(r,x) > u(T,x) > $(r,x) = Q(r,x) 

(7.23) 

V(T,x) G (0,00) X (0,00)". 

We have established the following result. 

Proposition 4. Recall the functions u(-,-), <^'(-,-) and Q{-,-) defined 
on (0,00) X (0,00)" by (4-1), (4-4) ^^^^^ C^-^V' respectively, and impose As- 
sumption C. Then (7.23) holds for every function U{-,-) £U. 

Remark 5. Here is a situation where Assumption C prevails: Suppose 
that (7.4) holds and that, for every z G (0, cxd)" and a G A{z), we have (9, a) G 
/C(z) for 9 given by 9iy = sj^, u = 1, . . . ,n and ss' = a. Then for any 

progressively measurable a:[0,oo) x — )• S*^ that satisfies (7.21) we select 
the progressively measurable functional [0,cxd) x — )• M" via 'di,{t,Lo) = 
Yl^=i CTjuit, (jj), u = 1, . . . ,n. This choice induces 

Mt,u^) = Ey- ^,(t) + . +^„,(t) J^-(^'^)' z. = l,...,n 

which obeys \\d{t, tr) f dt < 00 as in (7.22) for all (T, tr) G [0, 00) x SB; the 
process A(-) of (7.17) and (7.18) is a Q-martingale, whereas (2.4) becomes 



dXi{t)=X,{t) 



J2<^Ut,^)^W^it) + ^aij{t,X) dt 

M=l \j=l / . 



Xi{0)=Xi>0 



for i = 1, . . . ,n. The condition (7.21) guarantees now that X(-) takes values 
in (0,00)", P-^-a.s. in the resulting primal admissible system A4 G 9Jt(x). 
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8. Dynamic programming. The quantity Q{T,x.) defined in (7.11) is the 
value of a stochastic control problem: namely, the maximal containment^^ 
probability, over all measures with J\f S OT(x), that the process X{-) 
with dynamics (7.7), initial configuration X(0) = x G (0,00)*^, and controlled 
through the choice of progressively measurable functional a(-, •) as in (7.8), 
(7.10), does not hit the boundary of the positive orthant by time T. 

Let us suppose that the resulting function Q(-, •) is continuous on (0, 00) x 
(0,00)". Then it can be checked [as in Lions (1984), Lemma II. 1 and Lions 
(1983a), Theorem II.4] that it satisfies as well the following dynamic pro- 
gramming principle: for every initial configuration x G (0,00)", the process 

Q{T-t,X{t))l{r>t}, 

(8.1) 

<t<T is a -supermartingale, MM S ^R(x). 

[See El Karoui, Huii Nguyen and Jeanblanc-Picque (1987), Haussmann and 
Lepeltier (1990), Fleming and Soner (1993) and Krylov (1980) for results in 
a similar vein.] Equivalently, the process L{t)X(t)Q{T — t,X{t)), <t <T 
is a P-^-supermartingale, for every Ai G S[R(x). 

Consider now an arbitrary continuous function IJ : [0, 00) x (0, 00)" — t- 
[0,00) which satisfies U{0,-) = 1 on (0,oo)'", and is such that, for every 
Af G OT(x) and x G (0,cx))", the process U{T - t,X{t))l{r>t},0 <t<T is 
a Q^-supermartingale. We shall denote by ^ the collection of all such func- 
tions and note that C ^ and Q G From optional sampling we have then 
for every M G 0^(x) the comparisons 

U{T, x) > E< [^7(0, X(r))l|r>T}] = Q^iT > T), 

thus also c/(r,x) > (5(r,x), V(r,x) G (0,00) X (0,00)". 

In other words, the function Q{-, ■) defined in (7.11) is the smallest element 
of the collection ^. It is also clear from this line of reasoning that Mo G ^(x) 
attains the supremum Q{T,x) = sup_^gf)^(x) Q^(7' > T) in (7.11), if and only 
if the process Q{T - t, X(t))l{7->t}, < t < T is a -martingale. 

Theorem 1. Suppose that Assumption C and conditions (2.9)-(2.11) 
hold and that the function Q{-, •) of (7.11) is continuous on (0, 00) x (0, 00)". 
Then the infimum in (4-1) is attained, and 

(8.2) u(r,x) = $(r,x) = g(r,x) v(r,x) g (o,oo) x (0,00)^^. 

Proof. Consider an arbitrary function U{-,-) in the collection ^ just 
defined and fix an arbitrary pair (T, x) G (0, oo) x (0, oo)"; then for every e > 
0, consider a mollification Ue(-,-) GZ^ of the function U{-,-) with Us(T,x) < 
C/(T,x) + e. 
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Proposition 2 gives then u(r, x) < [/(T, x) +e. Since e > is arbitrary, 
this shows that u(T, x) is dominated by Q(T, x), the infimum of U (T, x) over 
aU functions U{-,-) ^ But the reverse inequality u(T, x) > Q(T, x) holds 
on the strength of (7.23), so (8.2) follows. □ 

8.1. The HJB equation. Under the conditions of Theorem 1, the arbi- 
trage function u(-, •) is equal to the function (5(-, •) of (7.11) and is continuous 
on (0,oo) X (0,00)". Thanks to the dynamic programming principle of (8.1), 
it is also a viscosity solution of the Hamilton-Jacobi-Bellman (HJB) equa- 
tion 

/ N dU , , v^v^ A ^2../ N DiU(T,z) \ 

(8.3) ^(r,z)= sup ^^z,z,a,J-DlUiT,z)+ ^ ) 

on (0,00) X (0,00)" [cf. Lions (1984), Theorem III.l or Lions (1983b), The- 
orem I.l]. 

If in addition to being continuous, as we assumed in Theorem 1, the 
function Q{-,-) of (7.11) is of class C^'^ locally on (0,oo) x (0,oo)", then 
the arbitrage function u(-,-) is not only a viscosity solution but actually 
a classical solution of the HJB equation (8.3). This is the case, for instance, 
under the combined conditions of Theorem 1 and Proposition 3; then the 
arbitrage function u(-,-) also satisfies on the domain (0,oo) x (0,oo)" the 
linear parabolic equation 

(8-4) fir,.) = ±±z.z,.,iz)(lDf^Uir,z) + ^^^) 
i=l j=l ^ 1 -r -r n/ 

with a: (0,oo)" — )• as in Assumption A or Remark 2, in addition to the 
initial condition 

(8.5) U{0,-) = 1 on (0,00)''. 

In particular, the arbitrage function u(-, •) satisfies, in this case, the require- 
ment (6.10) and belongs to the class U of Section 5. 

Recalling Propositions 1-3 and Theorem 1, we summarize the above dis- 
cussion as follows. 

Theorem 2. Suppose that conditions (2.9)-(2.11), (6.1), (6.2) and As- 
sumptions A, B and C are in force. 

Then the arbitrage function u(-, •) is the smallest element of the class U, 
as well as a classical solution of both the HJB equation (8.3) and of the 
linear parabolic equation (8.4), subject to (8.5). Furthermore (8.2) holds, 
the infimum in (4-1) is attained, and the Markovian investment rule tt^ {■, •) 
in (5.10) with U = \x satisfies (5.11) for every admissible system M G 9Jt(x). 
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Remark 6. We note that Theorem 2 is in agreement with general reg- 
ularity theory for fully nonlinear parabolic equations, as in Lions (1983c), 
Theorem II. 4 (see also Krylov (1987), Section 6.5; Krylov (1990); Wang 
(1992a, 1992b, 1992c), Theorems II.3.2 and III.2; or Lieberman (1996), 
Chapter XIV). 

As we mentioned already. Assumptions A and B can be replaced in Theo- 
rem 2 by the conditions of Remark 2. We conjecture that the conclusions of 
Theorem 2 should hold under even weaker assumptions but leave this issue 
for future research. 

We also remark that the function V{t,z) := (zi + • • • + z„)u(t,z), {t,z) £ 
(0,oo) X (0,oo)" satisfies an HJB-type equation simpler than (8.3), namely, 
the Pucci maximal equation, 

dV 1 n n 

(8.6) — (r,z) = - sup V V2;i2ja»j-D^jF(T,z), 

^ aeA(z) 

along with the initial condition 1/(0, z) = zi + ■ ■ ■ + Zn- In the setting of 
Theorem 2, equation (8.6) reduces to 

^(^,z) = -^^ZiZjaij(z)i?2y(r,z). 

i=l 3=1 

8.2. An example. Let us go back to the volatility-stabilized model in- 
troduced in Fernholz and Karatzas (2005), but now with some "Knightian" 
uncertainty regarding its volatility structure 

l<aii{t)^ii{t)<l + 5, 0<t<oo 

for some given (5 > 0. The case 5 = corresponds to the variance structure 
of the model studied in Fernholz and Karatzas (2005, 2009). 

More specifically let us assume that, for any given y G (3„, the compact, 
convex subset ^(y) of in (1.6) consists of all matrices a = {aij}i<ij<n, 
with aij = for j / i and 

yiO-ii = ri^{yi H hy„); i = 1, . . . ,n, 1 < r? < 1 + 5. 

The sets of (1.1) are given as 

/C(y) = {(0, a)|a G ^(y), ^ = (C^^, . . . , CV^)' with C G [v^, v^]} 

for some given constants Ci G (0, 1], C2 G (l,oo); these choices satisfy (7.4), 
(7.5). 

Condition (6.2) is satisfied in this case automatically (in fact, with A = 1), 
as are (6.1) and (6.6): it suffices to take 



(8.7) 



aii(z) = (2;iH Vzn)/zi, i = l,...,n 
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and -ff(z) = ^1,1=1 log -2^1) which induces Oi{z) = Y^ajj(z) in (6.3). These func- 
tions are ah locally bounded and locally Lipschitz continuous on (0, oo)". 
The HJB equation (8.3) satisfied by the arbitrage function u(-, •) becomes 



du 

— (r,z)= sup 
OT l<r)<l+<5 



n n 

- ^^(^1 + • • ■ + Zn)ziDlU (r, z) + Y,ZiDiU (r, z) 

i=l i=l 



and reduces to the linear parabolic equation 

du 1 ^ ^ 

(8.8) — (r, z) = - Y^i^i + ■■■ + Zn)ziDlU{T, z) + ZiD.iU{T, z) 

i=l i=l 

of (8.4) for the choice of variances in (8.7). The reason for this reduction 
is that the expression on the left-hand side of (8.8) is negative, so we have 
considerable simplification in this case. 

Remark 7. In this example, the arbitrage function u(-, •) can be repre- 
sented as 

in terms of the components of the (0, oo)"-valued capitalization process 
X(-) = (Xi(-), . . . ,X„(-))'. These are now time-changed versions Xi{-) = 
^j(A(-)), i = 1, . . . ,n of the independent squared-Bessel processes 



d^,{u)=Audu + 2^/W~(^dPi{u), mi{{)) = zi, 

run with a time change A{t) = (1/4) J^{Xi{s) + ••• + X„(s))ds common 
for all components, and with /3i(-), . . . , /3ri,(") independent standard Brown- 
ian motions [see Fernholz and Karatzas (2005, 2009), Goia (2009) and Pall 
(2011) for more details]. 

9. A stochastic game. For any given investment rule 11 € ^ and admis- 
sible system G 9K(x), let us consider the quantity 

(9.1) en,>((T,x) :=inf{r > 0:P^(Z^-^W'n(r) > X{T)) = 1}. 

This measures, as a proportion of the initial total market capitalization, the 
smallest initial capital that an investor who uses the rule 11 and operates 
within the market model M , needs to set aside at time t = in order for 
his wealth to be able to "catch up with the market portfolio" by time t = T, 
with P-^ -probability one. 

Our next result exhibits the arbitrage function u(-, •) of (4.1) as the min- 
max value of a zero-sum stochastic game between two players: the investor, 
who tries to select the rule 11 G so as to make the quantity of (9.1) as 
small as possible and "nature," or the goddess Tyche herself, who tries to 
thwart him by choosing the admissible system or "model" Ad G 9Jl(x) to his 
detriment. 
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Theorem 3. Under the conditions of Theorem 1, we have 



(9.2) u(r,x) = inf ( sup ^n,x(7',x))= sup ( inf ^n,A4(r,x; 

Proof. For the quantities of (9.1) and (4.7) we claim 

(9.3) en,A4(r,x)>UA4(r,x) V(n,7W)GqJx9Jt(x). 

Indeed, if the set on the right-hand side of (9.1) is empty, we have 
^ll,M (T, x) = oo and nothing to prove; if, on the other hand, this set is 
not empty, then for any of its elements r > the process IS 
a P-^-supermartingale, and therefore (4.6), that is, r > um{T,:x.), stih holds 
and (9.3) follows again. 

Taking the infimum with respect to 11 S on the left-hand side of (9.3), 
then the supremum of both sides with respect to G 9Jt(x), we obtain 



G(r,x):= sup infen,A^(T,x) > sup UA^(r,x) = $(T,x) 
xean(x) ^ne^p / m€W{x) 

from (4.7). The quantity G{T,x.) is the lower value of the stochastic game 
under consideration. 

In order to complete the proof of (9.2) it suffices, on the strength of 
Theorem 1, to show that the upper value 



G(r,x):=inf sup en,A4(r,x) >G(r,x) 
of this game satisfies 

(9.4) G(r,x)<u(r,x). 

To see this, we introduce for each given investment rule 11 G *p the quantity 
fln(r, x) := inf{r > : p-^(^rX(o),n^^) > X(r)) = 1, V7W G 93t(x)}; 

that is, the smallest proportion r > of the initial market capitalization 
that allows an investor using the rule 11 to be able to "catch up with the 
market portfolio" by time t = T with P-^ -probability one, no matter which 
admissible system (model) M might materialize. We have clearly 

(9.5) f)n(r,x)>u(r,x)ven,A^(r,x) V(n,7W)G*px9Jl(x), 
which leads to 

(9.6) u(r,x) = inf f)n(r,x) > inf ( sup en,A4(T,x)) = G(r,x) 
and proves (9.4). □ 
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9.1. A least favorable model and the investor's best response. Let us place 
ourselves now in the context of Theorem 2 and observe that Proposition 1, 
along with Proposition 2 and its Corollary, yields 

(9.7) en,A^,(r,x)>UA^„(r,x) = $(r,x) = eno,A4o(T,x) vn G ^, 

by virtue of (9.3) for M = Mo and of (5.11) for U = Here Mo is the 
"least favorable admissible system" that attains the supremum over 5[R(x) 
in (4.4), and XIq = vr* denotes the investment rule of (5.10) with [/ = <!>. 

In this setting, the investment rule IIo G attains the infimum 
infnetp f)n(7'>x) =u(T,x) in (9.6), and we obtain then 

fln„(T,x) =u(r,x) = $(r,x) =^n„,A^o(r.x) > en„,A^(r,x) 

(9.8) 

yM G wi(x) 

on the strength of (9.5). Putting (9.7) and (9.8) together we deduce 
^n,Xo(r,x) >u(r,x) =^n„,Xo(^'X) >Cn„,>!(7',x) 

y{U,M)e^ix 9Jt(x), 

the saddle property of the pair {Ilo,Mo) G *P x 9K(x). 

In particular, the investment rule Ho = vr* of (5.10) with [/ = $ is seen 
to be the investor's best response to the least favorable admissible system 
Mo G 9Jt(x) of Proposition 1, and vice-versa. In this sense the investor, 
once he has figured out a least favorable admissible system Mo, can allow 
himself the luxury to "forget" about model uncertainty and concentrate on 
finding an investment rule Ho G ^ that satisfies ^n,Xo(^, x) > ^UoMoC^t^): 
Vn G ^ as in (9.7), that is, on outperforming the market portfolio with the 
least initial capital within the context of the least favorable model Mo- 
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